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Abstract 

In this paper, we have explained transverse electro- 
optic effect by quantum-photonic model (QPM). This 
model interpret this effect by photon-electron inter- 
action in attosecond regime. We simulate applied 
electric field on molecule and crystal by Monte-Carlo 
method in time domain when a light beam is prop- 
agated through the waveguide. We show how the 
waveguide response to an optical signal with different 
wavelengths when a transverse electric field applied 
to the waveguide. 

Keywords: Quantum-photonic model(QPM), 

photon-electron interaction, attosecond regime, 
electro-optic (EO) effect, NPP. 



1 Introduction 

Electro-optic polymers are particularly interesting for 
new device design and high-speed operation [lj-[7j. 
Organic optical materials like MNA, NPP, MAP have 
a high figure of merit in optical properties in compar- 
ison with inorganic optical materials such as BBO, 
LiNbOa, [8J. 2-methyl-4-nitroaniline (MNA) and N- 
(4-nitrophenyl)-L-prolinol (NPP) have the highest 
figure of merit between organic nonlinear optical ma- 
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terials [9]. Thus they are used for electro-optic and 
nonlinear optic applications [8j- [2 Ij. The most effec- 
tive element in optical phenomena could be the re- 
fractive index. Because of the large electro-optic co- 
efficient of organic material, a certain amount of re- 
fractive index (RI) change can be realized with lower 
driving voltage than in other EO mater ials(l. 5 F//im 
for NPP,[10j). Ledoux et al. work on linear and 
nonlinear properties of NPP crystal. They first pro- 
posed a Sellmeier set for RI of NPP that is based 
on measured refractive indices [llj. Banfi, Datta and 
co-workers design a more accurate Sellmeier set for 
RI of NPP for nonlinear optical studies [13j-[15j. 
Their data is based on classic and macroscopic mea- 
surements and approaches. Some physicians explain 
RI in molecular bases for typical material [22] -[25]. 
Some authors have calculated RI of real material 
by quantum mechanical approach, [26]- [30] . Recently 
it is done some calculations and simulations on RI 
of liquid crystal by some softwares based on quan- 
tum mechanics (QM) [31]. In this paper, EO effect 
has been analyzed using QPM. We use QPM to ex- 
plain linear optical phenomena in molecular scales. 
This model gives us a constitutive vision about phe- 
nomena and real materials [32] -[35]. This approach 
is based on four elements: 1-Quasi-classic principle 
for justification of optical phenomenon in molecular 
scales, 2- Knowledge of crystal network and its space 
shape, 3- Short range intramolecular and intermolec- 
ular forces, 4- Monte-Carlo time domain simulation. 
In this model we suppose a laser beam is a flow of pho- 
tons while passing through single crystal film, inter- 



Figure 2: The crystal packing of NPP. 



Figure 1 : The molecular compound of NPP. 

acts with delocalization 7r-electron system of organic 
molecule and delays the photon in every layer [35] , 
By precise calculation of these retardation in every 
layer, we obtain RI in specific wavelength and explain 
EO effect too. We show that the phase retardation 
of input light with different wavelengths is distinc- 
tive, when it travels through the waveguide. The re- 
sults obtained from this method are well agreed with 
experimental data. Our favorite organic molecule is 
NPP. 

2 Crystal and Molecular Structure of 
our Favorite Organic Compound 

Organic molecular units and conjugated polymer 
chains possessing 7r-electron systems usually form 
as centrosymmetric structures and thus, in the 
electric dipole approximation, do not show any 
linear electro-optic and nonlinear second order 
optical properties. The necessary acentric may be 
provided by first distorting the 7r-electron system by 
interaction with strong electron donor and acceptor 
groups [^. In NPP molecule, nitro group acts 
as an acceptor and the other main groups on the 
other side of benzene ring acts as a week donor (see 
fig. 1). NPP (Ciii^i4A^203) (Fig.l) crystallizes 
in the solid state in an acentric monoclinic (with 
space group P2i) structure and their parameters 
are:a=5.261A°,b=14.908A°,c=7.185A° ,/3=105.18° 
and in the wavelength range of 0.5 to 2/im is 
transparent. The most interesting property of 
NPP crystal is the proximity of the mean plane of 




Figure 4: The bond lengths (in angstrom) and angles of NPP 
molecule. 



molecule with the crystallographic plane (101); the 
angle between both of these planes being 11°. Nitro 
group of one molecule in downward connects to 
Prolinol group in upper by hydrogen bonding. The 
angle between b orientation of crystal and N(l)-N(2) 
axis (charge-transfer axis) is equal to 58.6° [8j. Fig. 
2 shows the crystal packing of the NPP. For accurate 
and valid simulation, these properties and angles 
have to be exerted. For benzene molecule, benzene 
ring is a circle (see figl3]); Fig HI demonstrates the 
bond lengths and angles of the NPP molecule. As 
we see in this figure the bond lengths in benzene ring 
are not same. In our simulation, for similarity we 
consider an ellipse correspond to circle for electron 
cloud. We obtained £=0.26 for the ellipse of NPP 
from simulation, fig|5] shows the comparison of a 
circle and an ellipse with £=0.26. 




(a) 




(c) 




Superposition of (a), (b) and (c) 

Figure 3: Electron cloud for Benzene molecule that obtained from Huckle theory [37] and its approximation. 




3 Our Model for Electro-Optic Effect 

For a biaxial crystal that Ux Uy ^ riz the equation 
of index ellipsoid is 
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(1) 



Assume that single crystal film lies in x-y plane and 
light propagation is in the direction; therefore in the 
presence of an electric field the equation of index el- 
lipsoid by assuming crystal symmetry will become: 

(2) 

for NPP that Ex and Ey are transverse electric field 
components, [lOj. For MNA crystal Eq.(2) will be- 
come: 



(3) 



where Ex is only transverse electric field compo- 
nent [21j. With appropriate rotational transforma- 
Figure 5: The comparison between circles for Benzene tion, these relations can been simplified. In NPP, 

molecule electron cloud and NPP molecule electron cloud, (ap- and r22 and in MNA, Tn and r2i is large COef- 

proximately). As we see the electron cloud of ^-electron in ^^-^^^^ Therefore and Uy for NPP simplified to 

Tlx - {l/2)nlri2Ey and Uy - {l/2)n^r22Ey and for 
MNA Ux - {l/2)nlriiEx and Uy - {l/2)nlr2iEx re- 
spectively. The phase retardation F, with an applied 
electric field in a typical linear transverse electrooptic 
modulator will be obtained as follows, 38 : 



NPP has distorted a little. 



F = 



{nlr22-nlri2)Ey] (4) 
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for NPP and 

r = (t)y-4>x = ^■[ny-nx-^{nlrii-nyr2i)Ej;] (5) 

for MNA crystal. The total phase difference between 
two perpendicular polarization of light (in our exam- 
ple Ex and Ey), is 

A = Ao + (5 (6) 

where Aq is due to linear birefringence and 6 is due 
to linear electro-optic effect. In this case 5 is much 
smaller than Aq. 

In sub-micron space scales and sub-femtoseconds 
time scales, the optical constants loses its stabiliza- 
tion and classical equations for linear and nonlinear 
optical phenomena are not useful [45j. Now, we sug- 
gest a microscopic model for linear electro-optic phe- 
nomenon. If we have a monochromatic laser beam 
with frequency u and intensity I, then we may attain 
average photon flux from relation ^39j : 

now if we assume thin single crystal film of NPP in 
"6" direction of crystal (or "z" axis) radiated by a 
He-Ne laser with: A=633 nm, average power=10mw 
and beamwidth=20 microns, then from (|4]) average 
photon flux is equal to lO'^'^ photons /{s — cm?) that 
signifies in every second 10^^ photons arrive to each 
centimeter square. Moreover from data of crystal 
in subsection {A) in every 36.5 (A° ) on z direction, 
one NPP molecule exists. Therefore in every second 
36.5 X 10^ photons interact with any molecule or in 
other words in every 27ns, (with assuming only lin- 
ear optic phenomenon exist and nonlinear optic phe- 
nomenon do not exist, approximately. Because laser 
watt is not much), one photon interacts with any 
NPP molecules. In each interaction between photon 
and electron in every layer of crystal, we suppose a 
delay time equal to (ith layer of the crystal). Total 
delay time for m layers in crystal region is equal to: 

m 

1=1 

Consequently required time for photon transmission 
in L length of crystal is equal to r, achieved from 
relation: 

L nL L , , 



where Cq is velocity of light in vacuum. By using 
this relation, we can relate macroscopic quantity n 
to microscopic quantity r^. In biaxial crystals, r and 
consequently n depends on polarization direction of 
incident light. Because dipole- field interaction con- 
clusion is different for any direction of molecule. If 
Tx would be a microscopic delay for interaction of 
x-polarization field with dipole (or charge transfer 
action) and Ty would be a microscopic delay for in- 
teraction of ^-polarization field with dipole then the 
final phase difference between these two fields (named 
phase retardation) will be: 

A(l) = iJ ■ {Tx - Ty) = U ■ ^{Tx, - Ty^) (9) 

Of course this relation give Aq of (6). We justify 5 
from our model in later subsections. 

4 Photon-Electron Interaction in At- 
tosecond Regime 

For analyzing the interaction of the electric field 
of the photons with 7r-electron system, the time- 
dependent Schroedinger equation has to be used: 

H{r,t).^{r,t)=ih^^^^ (10) 

with the Hamilton operator H representing the total 
energy of the matter-light system and the wave func- 
tion ^ representing the quantum state of this system 
with all detailed spatial and temporal information of 
all particles in it. First, the stationary Schrodinger 
for matter without any external interaction is usually 
applied: 

Hmatter{r).(pm{r) = Em'^ni{r) (11) 

The interaction of the photon field with 7r-electron 
system can be described by first-order perturbation 
theory. The Hamiltonian of (10) is split into the 
material steady-state Hamiltonian of (11) and the 
Hamiltonian of the interaction as a small disturbance: 

Hiv^t) = HYno^iiQrp -\- H inter actionif) (-^^) 

With this equation the temporal change of the co- 
efficient describing the transitions of the the parti- 
cle under the influence of the light can be calculated 
from: 
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with the integration over the whole volume V of the 
wavefunctions. The probability of the population of 
state p is given by the square of Cp and the transition 
probability uOp^rn for the transition from state m to 
state p is given by: 



d_ 

dt 



Cp{t)\^ (Xfi' 



2 



(14) 



wich is proportional to the square of the transition 
dipole moment iip^rn'- 



[ ^l-H, 
Jv 



interaction 



• ^mdV (15) 



The interaction operator is given for a one-electron 
system in the dipole approximation, assuming a ra- 
diation wavelength large compared to the dimension 
of the particle, by: 



H inter action {t^ — ^'^^ij' particle^ 



(16) 



with the electric charge e, the position of the particle 
center at Vparticie and r as the relative position of 
the charge from the particle center and the electric 
field vector E. For more general case, including large 
molecules the electric field can be better expressed 
with the vector potential A(r, t) which is source free: 

V •A{r,t) =0 (17) 

and the electric field follows from this potential by: 

1 d 



cdt 



A(r,t) 



and the magnetic field by: 

H{r, t) = rotA{r, t) 



(18) 



(19) 



With respect to the quantum description, the vector 
potential can be written as: 



A(r, t) = ^ e. 



h\r) 



(20) 



with the counter m for the different waves of light and 
thus of the electric field, as the direction of the 
field vector, as the wavelength of the light wave, 
V as the volume the waves are generated in and km 
as the wave vector of the mth wave. The bm and 
are photon absorbtion and emission operators which 
would be light amplitudes in the classical case. These 
operators fulfill the following relations: 

^m-^p ^p -^m ~ ^mp 



hm-hp - bp.bm = 



bt-b±=0 



(21) 



which result in the description of the energy of the 
electrical field by a sum over harmonic oscilators as: 



H field = ^ bl^bmhVn 



(22) 



and the Hamilton operator for a single electron in 
the potential of the cores V and the electric field A 
is given by: 



Hele 



1 



2me^e 



-\p-eeA{r,t)f + V{r,t) (23) 



with the mass rueiectron and charge ee of the electron 
and the pulse operator: 



(24) 



With these definitions the interaction operator for a 
one-electron system follows from: 



interaction 



'^electron 



-A{r,t).p 



+ 



-A\r,t) (25) 



for linear interactions the second term can be ne- 
glected. But the interaction has to be considered for 
all charges in the particle which are in molecules for 
all molecules and core charges. The resulting inter- 
action operator is given by: 



Hinteractionij^^ t) — ^ ^ [ 

p 



'^electron 



A{rp,t).pp\ 



^core^q^e 
■^r.ore 



(26) 



with the charge Zcore.q of the qth core, the coordinate 
Rq of this core and its momentum Pq. In the dipole 
approximation the interaction operator for such a sys- 
tem can be written as: 

HinteracUon{r, t) = ^[-^£^i™Ii£i^(rp, t)] (27) 

^ '^charge,m 

and thus the transition dipole moment in the dipole 
approximation follows as: 



Zcharge,m^e 
'^charge,m 



E{rp,t)]).ipmdV 
(28) 
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For a real material such as NPP or MNA this formula 
would be very complicated and take enormous calcu- 
lations. Therefore some approximation must be ap- 
plied. It can be shown that for absorbtion or emission 
of photons the material has to perform a transition 
between two eigenstates Em and Ep of the material 
and thus the photon energy Ephoton has to fulfill the 
resonance condition: 

Ephoton — hVpiioton — \-^p -^m | (^^) 

But for our linear phenomenon, the photon en- 
ergy is about 2eV (in A = 630nm). If electron 
would be in HOMO (Highest Occupied Molecular 
Orbital), this electron do not go to LUMO (Low- 
est Unoccupied Molecular Orbital) or exited state 
by interaction. This phenomena is named nonres- 
onant phenomenon, [9l |40l [41] (nonresonant phenom- 
ena is not exclusive for nonlinear optical phenomena) . 
Therefor electron after interaction, go to quasi states 
that their life times is very short, then this electron 
go back to primary state after very short time. The 
nonresonant lifetime is determined by the uncertainty 
principle and the energy mismatch between photon 
energy in second time and the input photon energy. 
We can assume that the characteristic response time 
of this process is the time required for the electron 
cloud to become distorted in response to an applied 
optical field. This response time can be estimated as 
the orbital period of the electron in its motion around 
the nucleus which is about r 10~^^s or 100as,[41j. 
We can estimate this characteristic response time ac- 
cording to (9) if, n = 3, L = 3/i. Consequently Sr is 
equal to 10~^^sec. Because in b direction of crystal 
in 3/im length, approximately 4024 molecules exist, 
therefore the average quantity of r: 



is in order of 10~^^5 or 1 as. The perturbation in this 
very short time can assume semiclassically. In lin- 
ear phenomenon in this short time, just one photon 
interacts with one molecule. Because NPP molecule 
has delocalization electrons, (or 7r-electron system), in 
benzene ring, that photon interacts with this elec- 
tron type, [32] and it is annihilated [38j. We call this 
photon, a successful photon, (that does not produce 
phonon). 

To obtain 7r-electron wavefunction for benzene 
molecule the Schrodinger equation may be solved. 
Since this is very complicated process, it cannot be 



done exactly, an approximated procedure known as 
Hiickle method must be employed. In this method, 
by using Hiickle Molecular-Orbital (HMO) calcula- 
tion, a wave function is formulated that is a linear 
combination of the atomic orbitals (LCAO) that have 
overlapped [37] (see Fig. 3): 

* = ^Ci$i (30) 

i 

where the refers to atomic orbitals of carbon atoms 
in the ring and the summation is over the six C atoms. 
The coefficients of the atomic orbitals are calculated 
self-consistently through Roothaan's equations to ob- 
tain the ^ and the corresponding the one electron 
energies e, [43] : 

V 

where the Fock matrix F^.^ and overlap integrals S^^ 
are given by: 

Fy.. = + ^PAa[(M^|Aa) - \{^\\vg)\ (32) 

Act 

and 

Syv = y" <i)y4vdT (33) 

Here the core Hamiltonian matrix liy^v^ density ma- 
trix P^^, and two-electron repulsion integrals are 
given by: 

= I ^^H'^'^^dT, (34) 

occupied 

Pyv — 2 ^ ^ CyiCyi^ (35) 

and 

= [ (/)^(1)(/),(1) — (/)A(2)(/),(2)dridr2, (36) 
J ri2 

and the molecular Hamiltonian by 

i i<j ^ A ^"^^ 

(37) 

where the sum on i{A) is over electrons (nucleii), and 
the Za is the net core charge. The \Ci\^ is the prob- 
ability of the TT-electron at zth atom. Thus: 

\Ci? + \C2? + |C3p + + + ICep = 1 
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Figure 6: Assumed 7r-electron orbit of NPP molecule that 
obtained from fig. 3 for Benzene molecule (approximately). 




t = t„ e = e, 



Figure 7: 7r-electron system approximation for NPP 
molecule, ^ = is in the positive direction of dipole (Proli- 
nol side) and ^ = tt is in the negative direction of dipole (Nitro 
side). 



In the case of Benzene molecule: 




as followed from the symmetry of the ring [44l [45] . 
But NPP and MNA molecules aren't such as Ben- 
zene molecule. NPP is polar molecule. Nitro {NO2) 
is more powerful electronegative compound than pro- 
linol and pulls 7r-electron system; consequently, the 
probability of finding 7r-electron system at various 
carbon atoms of main ring isn't the same and the 
probability of finding 7r-electrons near the Nitro 
group is greater than near the prolinol group. There- 
fore there is no symmetry for NPP and electron cloud 
is spindly or oblong, (similar to dom-bell) (Figl6|). We 
estimate this form of electron cloud by an ellipse that 
our calculations would be uncomplicated. We assume 
effective positive charge that is located in one of fo- 
cal points of ellipse. The quantity of this effective 
positive charge is determined by semiclassical argu- 
ments. For attaining probability of electron presence 
on an orbit (FigCj), we say, T time is required by ra- 
dial vector to sweep total n.u.v interior area of ellipse 
(u and V are semimajor and semiminor axis of ellipse 




1 2 3 4 5 6 

Teta in Radian 



Figure 8: PDF of presence probability of 7r-electron in dif- 
ferent energy states when it rotates around the assumed orbit 
in FiglXl as we see the PDF is correlated to 0. 

respectively), in t times, this radial vector sweeps: 

t 

TT.U.V.— 

T 

area of ellipse, (see Figl7|). If t is the time, that elec- 
tron sweeps radian of orbit then t is obtained from 
this relation [32j: 

t = — {2arctan(W- tan(-)) 77;^^} 

27r^ + £ 2 l + 6:.cos((9) ^ 

(38) 

Where s is ellipse eccentricity. By using this rela- 
tion, we attain the required time (t) for electron to 
traverse from to -\- dO and it is divided by total 
time T. By this approach, we can determine the PDF 
(Probability Density Function) approximately. The 
PDF in apogee (near the Nitro group), is maximum 
and in perigee (near the Prolinol or Methane) is min- 
imum. Therefore PDF is correlated to from (38) 
and seen in Fig|8l The r quantity is correlated to 
and consequently, r quantity is correlated to presence 
probability of 7r-electron system. 
The angle between Y vector and charge transfer ac- 
tion (N1-N2) is 58.6° and X and Z axis is perpendic- 
ular to Y (Fig|9]). We consider propagation along Z 
direction. We spot a photon interacts with 7r-electron 
of NPP in first layer. After interaction, this photon 
gives its energy to electron and is annihilated. Elec- 
tron absorbs energy and digresses in direction of pho- 
ton momentum. Electron with photon energy, may 
not be unbounded and after arriving to apogee of di- 
gression, it returns back to ground state, because the 
photon energy is equal to hu = 1.96ev (h is Planck's 
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Figure 9: X, Y and Z axis and NPP molecule in dielectric 
frame. 

constant and v is frequency of laser beam) whereas 
energy for excitation is greater than Zev. When elec- 
tron returns to ground state one photon is produced. 
The time coming up and down is r delay time. This 
photon after freedom goes to second layer in direc- 
tion of annihilated photon (nonce, we assume the po- 
larization doesn't change), in second layer this pho- 
ton interacts with another delocalization 7r-electron 
certainly, because the effective interaction range of 
photon is approximately equal to its wavelength and 
is very greater than the distance between molecules. 
This molecule is nearest to photon effective central. 
This action is repeated for each layer. The location 
of photon-electron interaction is significant in every 
molecule and it is effective on r quantity directly. We 
assume that interacting photon has circle polarization 
and electron subject to virtual positive charge center. 
The phase retardation between and Ey (9) can be 
obtain from this relation, [32]: 



= ""-y^^ ^[cos(^,) - .^^m-rl (39) 

that 

~ l + e.cos(i9i) 

where e is the elliptical eccentricity and u is the 
semimajor axis of the ellipse. By applying an ex- 
ternal transverse electric field to organic crystal (in 
the range of several volts per micron) the shape of 
TT-electron system will be deformed slightly and we 
would be expect some noticeable variations in micro- 
scopic delay parameters (rx^Ty) and phase retarda- 
tion; (see Fig. 10). Consequently, by step-like change 
of input voltage, optical signal is switched between 
output port of DOS. 



5 Variation Analysis of Phase Re- 
tardation versus Applied Electric 
Field 

We simulate phase retardation of 3/im-length NPP 
crystal by Monte-Carlo method, then we generate 
random number using MAT LAB program. This pro- 
gram produces PDF quantities was explained in be- 
fore subsection and relates each of them to every 
molecule. These values are indexing -zr-electron po- 
sitions in each layer, by assumption a reference point 
(see fig. 6). Additionally we have used a MAT LAB 
program for Monte-Carlo simulation. The inputs of 
this program are: 

1. The wavelength of incident optical beam in 
which we want to design DOSs; 

2. /i,m,g, /c = 4:^77,^0 that are Planck's constant, 
electron rest mass, elementary charge. Coulomb con- 
stant and speed of light respectively. 

3. Unit cell parameters of NPP crystal: a, b, c, P 
and its other parameters that have given in subsec- 
tion (A). 

4. L: crystal thickness that in our simulation it is 
Sjurn. 

And the outputs of MATLAB program are: phase 
retardation in each wavelength. 

System calibration is done semiclassically by ex- 
perimental refractive index data.t In this method that 
we obtain three refractive indexes with x-polarization 
in threea with e (eccentricity), u (semimajor axis of 
ellipse) and Z (equivalent positive charge) in a way 
that refractive indexes in three wavelengths are very 
close to experimental data. Then we would see that 
refractive index in other wavelengths and other po- 
larization with same e, u and Z will be achieved. 
Of course these values, e, u and Z would be close 
to experimental structure of crystal, for example u 
would be greater than and smaller than minimum 
and maximum sizes of six lengths of benzene hexago- 
nal respectively, or e would be small but greater than 
zero. In other hand these values must be logical. 
From this method in our simulation we have obtained 
e = 0.26, Z = 3.9, u = 1.4A^ that is very close to ex- 
perimental and structural data. 

Xu and co-workers 26j , have done some electro- 
optic experiment about single crystal film of NPP. 
They have obtained |n^ri2 — ^^^22! = SAOpm/V and 
ri2 = 65pm /V in an optical beam with 1064nm wave- 
length. They have studied phase retardation between 
Ex and Ey of optical beam as a function of angle 
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Ellipse Eccentricity (External Electric Field) 



Figure 10: Phase retardation between Ex and Ey of optical signal with different wavelength by applied external electric 
field. As we see; by applied electric field, the virtual ellipse eccentricity change (electron in electrical field). With eccentricity 
variation, the microscopic delay (rx,ry) change and consequently phase retardation is modified. 



between electric field and charge transfer action of 
NPP. They have concluded that the maximum phase 
retardation was observed for the field oriented along 
the charge-transfer axis which was parallel to the film 
surface. The electro-optic effect or phase retardation 
was negligible when the electric field was applied per- 
pendicular to the charge- transfer axis. This concept 
could be justified by our model in previous subsec- 
tion. When the angle between charge-transfer action 
and external electric field change, the ellipse eccen- 
tricity modify and consequently, the phase retarda- 
tion alter. Obviously, from Fig. 5 when the external 
electric is parallel to charge-transfer axis, the ellipse 
drag more and the ellipse convert to a line. Therefor, 
ellipse eccentricity arise and from Fig. 9 the phase re- 
tardation growth. In the other hand, from Fig. 5 when 
the angle between external electric field and charge- 
transfer axis change, the ellipse is gathered and con- 
vert to circle and the eccentricity decrease to zero. 
Thus from Fig. 9 the phase retardation is lowered. 



6 Conclusion 

we justified linear EO phenomenon by QPM. This 
suggested physical model could be a powerful tool 
for analyzing and explaining processes that happen 
in waveguides with microscopic and nanoscopic sizes. 
We showed how the phase retardation between dif- 
ferent arguments of an optical field with distinctive 
wavelengths can take place. 
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